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ABSTRACT 

Telescopes often modify the input polarization of a source so that the measured circular or linear 
output state of the optical signal can be signficantly different from the input. This mixing, or polar- 
ization "cross-talk", is defined by the optical system Mueller matrix. We describe here an efficient 
method for recovering the input polarization state of the light and the full 4x4 Mueller matrix of 
the telescope with an accuracy of a few percent without external masks or telescope hardware mod- 
ification. Observations of the bright, highly polarized daytime sky using the Haleakala 3.7m AEOS 
telescope and a coude spectropolarimeter demonstrate the technique. 

Subject headings: Astronomical Techniques, Astronomical Instrumentation, Data Analysis and Tech- 
niques 



o 



> 

in 

vn 
o 



1. INTRODUCTION 

Spectropolarimetry is a powerful tool which may be 
enhanced by efficient techniques for telescope-detector 
polarization calibration. Unfortunately, it is common 
for modern altitude-azimuth telescopes with multiple 
mirror reflections to scramble the input polarization 
of the light before it reaches the analyzer of a coude, 
Nasmyth, or Gregorian polarimeter. Solar astronomers 
in particular have developed a repretoire of tech- 
niques for c a librat ii ig telescope s (lElinore et al. 1 119921: 
Kuhn et al.1 Il993t iGiro et al. I [20031: iSocas-Navarro ' 
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llchimoto et al. I l2005 



lElmore et al. I 120101 O 



ur ex- 



Socas-Navarro et al. 
Keller fc Snik I I200c 
perience with the Advanced Electro- Optical System 
(AEOS) 3.67m telescope is that linearly polarized 
light can even, in some circumstances, be transformed 
to nearly 100% circularly polarized light. There are 
7 highly oblique reflections before polarization is 
analyzed at our High Resolution Visi ble and Infra- 
red Spectropolari meter ( Hi VIS) (cf. iHodapp et al. 



^jThornton et al. I 120031 



iHar i 



Har rington et al. 
I2007L 



2006', '2009, 2010, 'Harringto n fc Kuhn 
Harringtion & Kuhn 2009a. b). 

Effective calibration of the telescope and instrument 
minimizes these systematic errors and enhances the po- 
larimetric precision of an observation. Moving optical 
elements, detector effects such as non-uniform sensitiv- 
ity or nonlinearity, atmospheric seeing, and transparency 
variations can all induce polarization errors that are of- 
ten larger than the photon noise. Strategies to minimize, 
stabilize and calibrate the instrument and telescope of- 
ten involve use of non-moving optics and emphasize rapid 
chopping and common-path differential techniques. Cali- 
bration techniques include observing unpolarized and po- 
larized standard stars, and use of polarized optical cal- 
ibration units that are designed to mimic the nominal 
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illumination of the telescope optics during target obser- 
vations. 

The new technique presented here is designed to yield 
the full Mueller matrix calibration of many telescope- 
polarimeter systems with an accuracy of a few percent. 
It requires only two or more observations of the linearly 
polarized sky due to solar Rayleigh scattering. The day- 
time sky is normally easily observable, highly polarized, 
and is a relatively well characterized calibration source. 
It is observable without using precious dark observing 
time and requires no hardware calibration modifications 
of the telescope. As long as multiple linear polarization 
input states are detected at each telescope pointing, the 
full telescope Mueller matrix can be recovered with an ac- 
curacy of a few percent. This is realized by observing the 
solar illuminated sky with the Sun at different zenith and 
azimuth positions. This paper illustrates the technique 
using the Haleakala Air Force AEOS 3.7m telescope on 
Maui. 

1.1. Polarization 

The followin g discussion of polarization fo rmalism 
closely follows iCoUet I [199^ and iClarke I [2001 In the 
Stokes formalism, the polarization state of light is de- 
noted as a 4-vector: 



S, = [/,Q,{/,y]^ 



(1) 
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In this formalism, / represents the total intensity, Q 
and U the linearly polarized intensity along polarization 
position angles 0° and 45° in the plane perpendicular to 
the light beam, and V is the right-handed circularly po- 
larized intensity. Note that according to this definition, 
linear polarization along angles 90° and 135° will be de- 
noted as —Q and — J7, respectively. The typical conven- 
tion for astronomical polarimetry is for the -\-Q electric 
field vibration direction to be aligned North-South while 
-\-U has the electric field vibration direction aligned to 
North-East and South- West. 

The normalized Stokes parameters are denoted with 
lower case and are defined as: 
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(2) 



The degree of polarization can be defined as a ratio of 
polarized light to the total intensity of the beam: 



P 



v/Q2 + c/2 + y2 



(3) 



For details on p olarization o f lig ht and stellar spec- 
tropolarimetry, see lCoUet lll992l and lClark e" 2009. 

To describe how polarized light propagates through 
any optical system, the Mueller matrix is constructed 
which specifies how the incident polarization state is 
transferred to the output polarization state. The Mueller 
matrix is a 4 x 4 set of transfer coefficients which when 
multiplied by the input Stokes vector (Si-^^^j) gives the 
output Stokes vector (Si^^^p^J: 



(4) 

If the Mueller matrix for a system is known, then one 
inverts the matrix and deprojects a set of measurements 
to recover the inputs. One can represent the individ- 
ual Mueller matrix terms as describing how one incident 
state transfers to another. In this work we will use the 
notation: 



M,: 



/ // QI UI VI' 
IQ QQ UQ VQ 
lU QU UU VU 

\IV QV UV VV 



(5) 



1.2. 



Deriving Telescope Mueller Matrix Elements 

Typical calibration schemes on smaller telescopes of- 
ten use fixed polarizing filters placed over the telescope 
aperture to provide known input states that are de- 
tected and yield terms of the Mueller m atrix. This ap- 
proach has been used for solar telesco pes (|Socas- Navarro I 
I2005al|bl . iSocas-Navarro et al. I[2005ll2006l ). Another ap- 
proach, also used in the solar observations has been 
to image known sources and use spectropolarimet- 
ric data and Zeeman effect symmetry properties to 
isolate and determ i ne terms in t he M ueller matrix 
(jKuhn et al. I I1993L lElmore et al. 1 120101 ). Night-time 
observations can use unpolarized and polarized stan- 
dard stars to measure polarization properties of tele- 
scopes fcf. IHsu & BrcKcr 1982, Schmidt ct al. 1991 
iGil-Hutton fc Benavidez I 120031 IFossati et al. I 120071 ). 
Many studies have either measured and calibrated tele- 
scopes, measured mirror properties or attempted to 
design instrument s with min imal pola rimetric defects 
(cf. Sanchez Alm eida et aTl [T 991. Gi ro et al. I 120031 
Patat & Ronianiello ' '2006";Tinbcrgcn 20 071 iJoos et al. I 
200a. ,van Harte n ct al. 2009, Roelfscm a et al. II2O1O0 . 



1.3. The Polarized Sky as a Calibration 

The observed polarization of the daytime sky is a 
useful calibration source. Scattered sunlight is bright, 
highly (linearly) polarized and typically illuminates 
the telescope optics more realistically than calibration 
screens. A single-scattering atomic Rayleigh calcula- 
tion (cf . iCoulson ] 119801 11998D is often adequate to de- 
scribe the skylight polarization. More realistic model- 
ing, involving multiple scattering and aerosol scattering. 



is also available using industry standard atmospheric ra- 
diat ive transfer software packages such a s MODTRAN 
(cf. iFetrow et al~l [20021 [Berk et al. I [20061 ). 
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Fig. 1. — The celestial triangle representing the geometry for 
the sky polarization at any telescope pointing. 7 is the angu- 
lar distance between the telescope pointing and the sun. On 
is the solar angular distance from the zenith. 9 is the angular 
distance between the telescope pointing and the zenith. (j> is 
the angle between the zenith direction and the solar direction 
at the telescope pointing. 1/) is the angle between the tele- 
scope pointing and the solar direction at the zenith. The law 
of Cosines is used to solve for any angles needed to compute 
degree of polarization and position angle of polarization. 



There are many atmospheric and geometric con- 
siderations that determine the skylight polarization 
at a particular observatory site. The linear po- 
larization amplitude angle can depend on the so- 
lar elevation, atmospheric aerosol content, aerosol 
vertical distribution, aerosol scattering phase func- 
tion, wavelength of the obs ervation and seconda ry 
sources of illu mination (cf. iHorvath et al. I l2002al[b , 
Lee 1998, Liu fc Voss I [l997L ISuhai fc Horvath I I2004i 
Gal et al. 2001, 'Vcrmeulen ct al. "2000', 'Pomozi et al.J 
2001, Cronin et al. 2005, 2006, Hcgcdiis et al. 2007). 
Anisotropic scattered sunlight from reflections off land 
or water can be highly po larized a nd temporally 



variable ( Litvinov et al. I I2OIOI [P cltonic mi et al 



He et al. '2010', 'Salinas 
Kissclcv & Bulgarclli 
tical properties and 



2009 



2010 



fcf. IWufc.Tinj 



Vermeulen et al. 



1997 



L iew li2007nOta et al, 
2004j). Aerosol particle op- 
vertical distributions also vary 



iShukurov fc: Shukurov I 



2OOOI lOugolnikov fc Maslov I 



2006 



2002 



2005allbl . T2009ailbr lUeolnikov et al. II2004D . The polar 
ization can change across atmospheric absorption bands 
or can be influenced by other scat tering mechanisms 
(cf. IBoesche et 311120061 IZeng et al. , 2008 A ben et al. I 
1999, 200f). 

Deviations from a single scattering Rayleigh model 
grow as the aerosol, cloud, ground or sea-surface scat- 
tering sources affect the telescope line-of-sight. Clear, 
cloudless, low-aerosol conditions should yield high lin- 
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ear polarization amplitudes and small deviations in the 
polarization direction from a Rayleigh model. Obser- 
vations generally support this conclu sion (iPust fc Shawl 
[T005, 2006a^b, 2007, 2008, 2009, Sha w et al. \\2mW . 

The geometry of our Rayleigh sky model is seen in Fig- 
ure [T] The geometrical inputs are the observers location 
(latitude, longitude, elevation) and local time. The solar 
location and relevant angles from the telescope point- 
ing are computed from the spherical geometry in Fig- 
ure [TJ The maximum degree of polarization (Smax) in 
this model occurs at a scattering angle (7) of 90°. The 
Rayleigh sky model predicts the degree of polarization 
(5) at any telescope pointing as: 

= -, n 6) 

1 \- COS^^ ^ ' 

Since the angle of polarization in the Rayleigh sky 
model is always perpendicular to the scattering plane, 
one can derive the angle of polarization with respect to 
the altitude and azimuth axes of the telescope. The law 
of Cosines for the scattering plane breaks down at the 
Zenith where is 0° and one must simply calculate the 
difference in azimuth between the telescope pointing and 
the solar azimuth to find the orientation of the polariza- 
tion. An example of the model we compared with our 
observations is shown in Figure [5] for Haleakala on Jan- 
uary 27*'* 2010. 

One set of all-sky polarization measurements obtained 
at the Mauna Loa Observatory, a nearby site and at 
very similar altitude to Haleakala (jPahlberg et al. II2009D 
are particularly relevant. The maximum degree of po- 
larization {8raax) chaugcd throughout the day in their 
observations from a value of ^-^60% with the sun at 
an elevation of 30°, to 85% with the sun setting. Al- 
though this is a significant change in the maximum 
degree of polarization [bmax) presumably due to mul- 
tiple scattering, the polarization direction at all tele- 
scope pointin gs is expected to be cl o se to the Rayleigh 
mode l (cf. jPahlberg et al. I 120091 iSuhai fc Horvath I 
[2OOI IPiIsTfe Shaw 1120091 IPomozi et al. II2001D . In our 
analysis below we will use the Rayleigh model to define 
the input polarization angle. We believe our assumed 
input angle is a ccurate over most of the sky to <5° (eg. 
ISuhai fc Hof^th 2004 ). 

1.4. Using the Sky for Full-Stokes Calibration 

Despite the fact that the Rayleigh sky only provides 
linearly polarized light to the telescope, with reasonable 
assumptions, we can recover the full Mueller matrix. We 
find that the telescope and spectrograph only weakly po- 
larize unpolarized input light, and the optical system has 
only a weak depolarization effect on polarized input light. 
Thus, the first row and column of the system Mueller ma- 
trix can be approximated by the corresponding row and 
column of the identity matrix. Direct measurements of 
these terms using stellar observations with AEOS con- 
firm that over the visible spectrum these Mueller matrix 
terms are less than 5%. With this assumption the tele- 
scope system scrambles the input polarization simply by 
rotating the input 4 element Stokes vector in the 3-space 
defined along the Q, U , and V polarization axes to yield 
the output measured state. 




Fig. 2. — This Figure shows the sky polarization modeled for 
Haleakala at 3:00 UT on January 27*'' 2010. The degree of polar- 
ization ((5) at all telescope pointings with 5max scaled to 100% is 
in the top left panel. The scattering angle with respect to the sun 
is in the top right panel. The solar altitude is 15° and the solar az- 
imuth is 243° . The calculated q and u inputs as seen by the AEOS 
altitude-azimuth telescope are on the bottom two panels. -|-Q is 
defined as E-field vibrations along the + Altitude direction. +[/ is 
defined as E-field vibrations along the -|-Altitude, +Azimuth di- 
rection. All panels show stereographic projections to the observed 
sky with the Zenith in the center of the circle. North, East, South, 
West are on the perimeter of the circle with North up and West 
right. 

In general only two measurements with different in- 
put linear polarization orientations are sufficient to de- 
termine the telescope polarization properties. In prac- 
tice wc define a least-squares problem that takes several 
input polarization measurements of the sky at different 
times but with identical optical configurations to derive 
telescope properties. 

2. POLARIZED SKY OBSERVATIONS 

To demonstrate this technique, we collected sky obser- 
vations with AEOS using the new low spectral resolution 
mode for our coude spectrograph we call LoVIS. This 
new spectrograph has been characterized to verify the 
expected performance using standard stars, calibration 
optics and various tests outlined in the Appendix. The 
LoVIS system polarization properties most relevant to 
this work include the telescope and instrument induced 
polarization, depolarization, the polarimetric response of 
LoVIS. These are briefly described below and in the ap- 
pendix. 

The measured degree of polarization for the many Lo- 
VIS sky observations we have obtained shows that there 
is relatively low instrument depolarization. For exam- 
ple, on January 27*'' 2010 UT we have twenty sky po- 
larization observations with the telescope pointed at the 
Zenith taken over three hours ending just after sunset. 
The measured degree of polarization averaged over all 
wavele ngths is shown in Figu re [3] as the solid line. Fol- 
lowing iDahlberg et al. 1 120091 Fi gure 7, we assume the 
maximum degree of polarization (5max) for the Rayleigh 
sky as roughly 65% when the solar elevation is 30°, 75% 
at 15° and 85% at 0°. If we apply this simple linear rela- 
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Fig. 3. — This Figure shows the average sky polarization with 
the telescope pointed at the Zenith measured by LoVIS on January 
27*'' 2010 UT. The statistical error in the measurement is much 
smaller than the width of the line. The solid line is the average 
polarization measured for all wavelengths. The dashed line is a 
Rayleigh Sky model prediction for the polarization at the Zenith 
following measurements from Mauna Loa in Dahlbcrg et al. 2009 
Figure 7. The maximum degree of sky polarization in the 90° 
scattering plane {5max) is assumed to rise linearly from 65% to 
85% as the sun sets from elevation 30° to 0° . The scattering angle 
(7) between the sun and the Zenith telescope pointing is shown in 
the right hand y-axis. 

tion for 6max and compute the Rayleigh sky polarization 
(S) at the Zenith for our observations, we get the dashed 
curve in Figure [31 The Figure also shows the scattering 
angle (7) between the Zenith telescope pointing and the 
sun. 

The computed degree of polarization changes with 
wavelength for all twenty Jan 27*'' observations is shown 
in Figure [H The average polarization shown in Figure 
[3] was removed from the spectra leaving only residual 
chromatic changes. There is a mild trend with wave- 
length for the measured polarization to decrease by about 
10% from 55OOA to 75OOA. The atmospheric oxygen 
absorption band shows a mild increase in polarization. 
There is also a substantial change seen the chromatic 
trend seen in the final observation taken just after sun- 
set as shown by the dashed line of Figure |4] consistent 
with stu dies of the occultation and c hangin g illumin a- 
tion (cf. [O ugolnikov fc_Madovj[2QQl l2n?)5al lhl. [2009al lbl. 
lUgolniko7 et al. 2004). 

Given that AEOS telescope and LoVIS spectropo- 
larimeter induces polarization of less than 3% and the 
depolarization is of the same order, we can use a 3 x 3 
Mueller matrix approximation. In this 3x3 QUV space, 
the degree of polarization is preserved for any input state 
and the 3x3 Mueller matrix would be approximated as a 
rotation matrix inside the Poincare sphere. 

3. TELESCOPE POLARIZATION AS ROTATION 

To model the the 3x3 Mueller matrix as a rotation 
matrix, we will use Euler angles. We scaled all our mea- 
sured Stokes vectors to unit length in order to remove 
the residual effects from changes in the sky degree of 
polarization, telescope induced polarization and depo- 
larization; and to put our measurements on the Poincare 
sphere. We denote the 3 Euler angles as {a, (3, 7) and use 
a short-hand notation for sines and cosines where 003(7) 




5500 6000 6500 7000 
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Fig. 4. — This Figure shows the continuum polarization variation 
from average as a function of wavelength for all twenty LoVIS 
Zenith measurements on Jan. 27*'' 2010. The average value of 
the polarization, shown in Figure [3] has been removed from all 
observations to highlight the change with wavelength. The solid 
lines show observations taken 2-3 hours before sunset. The dot- 
dash lines show observations taken 0.2-1 hour before sunset. The 
dashed line shows a single observation taken just after sunset. 

is shortened to Cj. We specify the rotation matrix (My) 
using the ZXZ convention as: 

/ s.^ \ / 1 \ / Ca Sa \ 

My = -S7 C;3 S/3 ~Sa Ca 

V 1 J \ -sp cfs J \ 1/ 

With this definition for the rotation matrix, we solve 
for the Euler angles assuming a fully linearly polarized 
daytime sky as calibration input. If we denote the mea- 
sured Stokes parameters. Si, as {qrmUmjVm) and the in- 
put sky Stokes parameters, Rj, as (qr,Ur,0) then the 
3x3 QUV Mueller matrix elements at each wavelength 
are: 

/ qm\ / QQ UQ VQ\ / qr\ 

S,^[ u„, ]^ My R, ^ [ QU UU VU ] [ ur ] 
\v,nj \QV UV VV J \ Oj 

(8) 

This set of equations has 6 variables and only 5 known 
quantities. We have no V input to constrain the VQ, VU 
and VV terms. Nevertheless, two input and output vec- 
tors on the Poincare sphere are sufficient to fully specify 
the full 3x3 rotation matrix. Thus, we use the fact that 
the sky polarization changes orientation with time and 
take measurements at identical telescope pointings sep- 
arated by enough time for the solar sky illumination to 
change. This yields an over-constrained solvable prob- 
lem for all six linear polarization terms in the Mueller 
matrix. 

When using this rotation matrix approximation for the 
telescope Mueller matrix, the Rayleigh Sky input Stokes 
parameters multiply each term of the rotation matrix to 
give a system of equations for the three Euler anlges (a, 
(3, 7). This system of equations can be solved using a 
normal non-linear least-squares minimization by search- 
ing the (a, P, 7) space for minima in squared error. This 
direct solution of this set of equations using standard 
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Fig. 5. — This Figure shows the six Mueller matrix elements 
versus wavelength at a telescope pointing of 90° altitude, 225° az- 
imuth. These were derived from three sky measurements on Jan- 
uary using the two-step solution. The UV term grows 
with increasing wavelength while the QU term decreases with wave- 
length. The other Mueller matrix elements remain near or 1 and 
show smaller chromatic variation. 
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Fig. 6. — This Figure shows the variation in Euler angles with 
wavelength at a telescope pointing of 90° altitude, 225° azimuth. 
The value for each Euler angle at 6560A have been removed for clar- 
ity [a=176.5°, ^=81.6°, 7=98.8°]. The change with wavelength is 
dominated by /9 which corresponds to the changes in the QU and 
UV Mueller matrix elements from Figure [5] The statistical noise 
is smaller than the width of the line. 



minimization routines is subject to several ambiguities 
that affect convergence using standard minimization rou- 
tines. For example, Euler angles are symmetric under 
the exchange of (a, 7) with ( -180° +a, -/3, -180°-H7). 
Solutions are also identical with multiples of 2tt. As a 
convenient approach to this least-squares problem, we 
used a different set of equations in a two-step solution. 
We first perform simple least squares solution directly 
for the six linear-polarization Mueller matrix elements. 
In the second step we fit a rotation matrix to the six 
Mueller matrix elements at each wavelength. As an ex- 
ample. Figure [5] shows the Mueller matrix elements de- 
rived with the Jan. 27*'' 2010 observations. The details 
of our methods for deriving Euler angles and an example 
of how one could plan sky calibration observations are 
outlined in the Appendix. 

3.1. Euler Angle Solution Properties 

The Euler angles are well behaved smooth functions of 
both wavelength and telescope pointing. For instance, 
the Euler angle variations with wavelength calculated 
from January 27*'' 2010 observations at a telescope point- 
ing of at 90° altitude and 225° azimuth are shown in 
Figure [51 The change with wavelength is dominated by 
/3 while smaller variation is seen in a and 7. For clarity 
in the Figure, the Euler angle values at 6560A have been 
removed so all curves overlap. Similar curves are seen for 
other telescope pointings. 

The Euler angles are also smooth functions of telescope 
pointing. We have observations on multiple days taken 
on a grid of 8 azimuths and 3 different altitudes. Figure[7] 
shows the derived Euler angles at 6560A as a functions of 
azimuth. There are three different curves corresponding 
to the three different telescope altitudes. We find that 7 
effectively absorbs the changing azimuth while a strongly 
varies with altitude. Since there is substantial telescope 
induced rotation, these curves show significant deviations 
from a purely geometrical rotation. 

The derived Euler angles in Figure [6] are quite repeat- 
able when derived from calibration observations taken on 



different days. For instance, we have two sets of sky ob- 
servations taken December 10*'' and 11*'' 2009 each con- 
sisting of thee observation sets spaced roughly one hour 
apart. The Euler angles derived from on different days 
agree to within 0.5° giving an estimate of the system- 
atic error limits. The Mueller matrix elements derived 
with observations from different days agree to better than 
0.01. The changes in /3 with telescope azimuth shown in 
Figure [7] are well above the systematic noise of roughly 
0.5° and are repeatable even when using many different 
combinations of calibration data on both days. Similar 
trends for Euler angle variation with telescope pointing 
are seen at other wavelengths. 




100 200 300 
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Fig. 7. — This Figure shows the variation in derived Euler angles 
at 6560A as a function of telescope pointing. The top panel shows a 
for all telescope azimuths at elevations of 90° (solid), 75° (dashed) 
and 55° (dot-dashed). The middle panel shows /3 and the bottom 
panel shows 7 for all azimuths at the same elevations as a using 
the same line scheme. This Figure clearly shows that a changes 
the strongest with elevation while 7 shows changes with azimuth. 
There is not much change seen in /3. 
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3.2. Calibrating Spectropolarimetric Observations 

Once telescope calibrations have been derived, observa- 
tions subsequently taken with the instrument at the same 
telescope pointing can be de-rotated. The calibrated po- 
larization measurements are then oriented in a reference 
frame on the sky corrected for all geometric, telescope 
and instrument induced rotation. We have done many 
experiments where sky observations on some days are 
used to calibrate sky observations taken on other days 
to examine the stability and repeatability of the calibra- 
tions. 

As an example, Figure |8] shows six individual sky spec- 
tropolarimetric measurements taken on December 10*'' 
and 11*'' 2009 at different times of the afternoon. All 
observations were taken at a telescope pointing of az- 
imuth 021° and altitude 70°. The input linear polariza- 
tion orientation changes with time and there is significant 
chromatism and cross-talk observed in all Stokes param- 
eters. The three observations from each day are taken at 
similar times and the observed polarization spectra are 
similar. 
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Fig. 8. — This Figure shows the measured spectropolarimetry 
with LoVIS while the telescope is fixed at altitude 70° and azimuth 
021° scaled to 100% polarization. Three observations were taken 
December 10*'' 2009 (solid lines) and another three observations 
were taken at similar times on December 11*'' 2009 (dashed Unes). 
The statistical noise is smaller than the width of the line. 



In order to illustrate a general method of telescope 
calibration, imagine a scenario where we derive a set of 
calibrations using observations from one day and then 
use this telescope calibration to de-rotate observations 
taken on another day at an identical telescope pointing. 
We will treat the December 10*'' 2009 observations as a 
calibration set and derive Euler anlges for all telescope 
pointings and wavelengths. With these Euler angles in 
hand, we will then treat observations taken on December 
11*'' 2009 as our science targets. The de-rotated science 
observations will be compared against the Rayleigh sky 
prediction. This scenario is exactly what a typical ob- 
server would perform during routine operations. The 
difference between the de-rotated December 11*'' obser- 
vations and the expected sky input Stokes parameters 
are shown in Figure [HI The residual error is typically 
less than 0.01 for any Stokes parameter showing that the 
de-rotated observations match the expected sky model 
with very high accuracy. 
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Fig. 9. — This Figure shows the difference between the de-rotated 
measured Stokes parameters and the predicted sky input Stokes 
parameters at altitude 70° and azimuth 021°. The residual errors 
are less than 0.01 showing very low error between measured and 
expected Stokes parameters. 

We get consistent results for Euler angles using many 
different combinations of observations taken between Au- 
gust 2009 and March 2010. Over this time period there 
were changes to the optical configuration such as re- 
mounting of the retarders, switching gratings and chang- 
ing LoVIS / HiVIS modes as would be expected during 
normal observing operations. If one derives a residual ro- 
tation angle between the de-rotated measurements and 
the predicted sky, the residual rotation angles are al- 
most always less than 3° when using any calibration set 
against any other observation set. If one uses calibration 
data where there are no optical changes between calibra- 
tion and observation, this resulting angular error is much 
smaller, typically less than 0.5° rotation. This allows us 
to estimate systematic error from optical configuration 
changes as a few degrees residual rotation. When observ- 
ing as calibrated, the errors are under 1° QUV rotation. 

4. DISCUSSION 

We have demonstrated a simple method for deriving 
telescope polarization properties and Mueller matrix el- 
ements. By using observations of the bright, highly po- 
larized daytime sky and taking observations at multiple 
times using identical optical configurations one can ex- 
tract the polarization response of the telescope. 

This method can be easily implemented via least 
squares methods. A direct least-squares solution for 
Euler angles and a two-step solution for the linear- 
polarization sensitive Mueller matrix elements and a ro- 
tation matrix fit to these elements have been imple- 
mented. Examples of how to optimally prepare calibra- 
tion observations and solve for the telescope response 
properties are detailed in the Appendix. 

The AEOS telescope and LoVIS spectropolarimeter 
have over 20 refiections, many at high incidence angles 
before the polarization analyzer. However, the major po- 
larization cross-talk effect is a simple rotations in QUV 
space that is readily derived to an accuracy of a few per- 
cent from sky observations. Observations taken on mul- 
tiple days show we can de-rotate LoVIS sky observations 
to the predicted Rayleigh sky with at least 0.5° accuracy. 

The polarized sky is a calibration source that is bright. 
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easily observable, highly polarized, simple to model, well 
characterized and easy to verify with additional instru- 
mentation. Unpolarized and polarized standard stars 
have several disadvantages as calibration sources. Stars 
are not available at all telescope pointings and require us- 
ing night-time observing time for calibration. Telescopes 
with complex optical trains such as AEOS require cali- 
bration at many pointings making the unavailability of 
stars at all pointings a serious problem. Solar instru- 
ments are typically unable to observe stars at all but 
can easily utilize this sky-based technique. Polarized 
standard stars normally have low degrees of polarization 
which adds time or noise constraints to such calibration 
efforts. This technique uses sky-light to illuminate; the 
telescope and instrument optics like the target observa- 
tions. Night-time observations benefit from the fact that 
day-time observations are used to calibrate, rather than 
using precious night telescope time for calibration. 
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5. APPENDIX: SOLUTIONS FOR EULER ANGLES 

There are a number of methods available to solve for 
Eulcr Angles and Mueller matrix elements given a set 
of sky polarization observations. Suppose we have mea- 
surements at different times but at identical telescope 
pointings. Multiplying out terms in equation [7] for the 
ZXZ-convention we get: 



■Q '-'7 



— SaSj + CaCfjCj SpC^ (9) 



-CaSp 



C/3 



Then equating Mueller matrix elements to rotation 
matrix elements could write the system of equations for 
the three Euler angles as: 



/ <3rl(+CQC^ — CfjSaS-,) + Mr 1 ( + C-, + Cq S-, ) \ 

qrli + SaSfi) + Url{-CaSf,) 

qr2{+CaC^ — C/jSaS^) + «r2 (+C-, + Cq ) 

IJr2( — C,3C_,Sq — CaS-,) + «,-2 (+Ca C.^ — S.^ ) 

qT3{+CaC~, — CfiSaS^l) + «r3 (+Ct So + Cq ) 

1 qrii — CfjC^Sa — CaS-,) + «r3 (+Ca — ) 



(10) 

This system of equations can be solved using a nor- 
mal non-linear least-squares minimization by searching 
the (a, /3, 7) space for minima in squared error. With 
the measured Stokes vector (S^), the Rayleigh sky input 
vector (Ri) and a rotation matrix (Mij) we define the 
error (e) as: 



fit a rotation matrix, we avoid over-polarizing. 

Provided the arrays are indexed properly, the normal 
solution for Mueller matrix elements can be computed 
via the normal least-squares method. We re-arrange the 
time- varying Rayleigh sky inputs to (Rij) for i indepen- 
dent observations and j input Stokes parameters. The 
measured Stokes parameters (Si) become individual col- 
umn vectors. The unknown Mueller matrix elements are 
also arranged as a column vector by output Stokes pa- 
rameter (Mj ) . If we write measured Stokes parameters 
as (jirai , WjTii , Vrai ) and the Rayleigh input Stokes parame- 
ters as ((jri , Uj. . ) , we can explicitly write a set of equations 
for just two Mueller matrix elements: 



S. = 




UQ 



(12) 



We have three such equations for each set of Mueller 
matrix elements. With this can express the residual er- 
ror {ti) for each incident Stokes parameter (S^) with an 
implied sum over j as: 



£4 



(13) 



The normal solution of an over specified system of 
equations is easily derived in a least-squares sense us- 
ing matrix notation. The total error E as the sum of all 
residuals for m independent observations we get: 



E 



(14) 



i=i j=i 

For n measurements, this gives us 3 x n terms in this 
least squares problem. The Euler angles give identical 
rotation matrix elements under the exchange of (a, /3, 7) 
with ( -180°-|-a, -180°+7) as well as with additional 
multiples of 27r. This direct solution is easily solvable 
in principle but the ambiguities make implementing this 
solution with existing software languages more cumber- 
some than necessary. 

As an alternative method to the direct least-squares so- 
lution for Euler angles, we can do a two-step process that 
gives the same result but is much easier to implement. 
First we solve a system of equations for the Mueller ma- 
trix elements directly that are not subject to rotational 
ambiguity. With the Mueller matrix elements in hand, 
we can then perform a simple fit of a rotation matrix 
elements to the derived Mueller matrix terms. This two- 
step process allows us to have accurate starting guesses 
to speed up the minimization process, resolve Euler angle 
ambiguities and to estimate error propagation. 

When deriving the Mueller matrix elements of the 
telescope, one must take care that the actual derived 
matrices are physical. For instance, there are various 
matrix properties and quantities one can derive to test 
the physicality of the matirx (iGiv ens & Kostinski 1993, 
iKostinski et al. 1 119931 iTakakura fc StoU ■ .2009 ) . Noise 
and systematic uncertainty might give over-polarizing or 
unphysical Mueller matrices such as the element above 
that is >1. By simply using Mueller matrix elements to 



We solve the least-quares system for the unknown 
Mueller matrix element (Mj) by minimizing the error 
with respect to each equation. The partial derivative of 
Equation [13] for with respect to Mj is just the sky 
input elements Rij. Taking the partial with respect to 
each input Stokes parameter we get: 



dE 



2 E = -2 E I^'^ - ^ R,.M. 



(15) 



We have inserted a dummy sum over the index fc. Mul- 
tiplying out the terms and rearranging gives us the nor- 
mal equations: 



Ry RifeMfe — ^ RjjSi 



(16) 



Which written in matrix notation is the familiar solu- 
tion of of a system of equations via the normal method: 



M 



R^S 
R^R 



(17) 



This simple equation is very easy to implement with 
a few lines of code provided your observation times are 
chosen to give you a range of input states for a well- 
conditioned inversion. The noise properties and inver- 
sion characteristics of this equation can be calculated in 
advance of observations and optimized. We can write 
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the matrix A with an implied sum over i observations 
for each term as: 



A 



R^R 



(18) 



The solution to the equations for the three sets of 
Mueller matrix elements can be written as: 



M, = 



QQ 
UQ 



M, 



( QU 
\ uu 



QV 

uv 



= A 



= A 



(19) 
(20) 
(21) 



As an example, if we compute the inverse of A and 
multiply out A^^ for the QQ term we can write: 



(grigm.)(Mri'»ri) - iUriqmi){qriUri) 



(22) 



In this manner, we can easily implement the usual ma- 
trix formalism with a time-series of daytime sky obser- 
vations to measure six Mueller matrix elements. 

With the Mueller matrix terms in hand, solving for 
accurate initial Euler angle guesses is straightforward. 
As an example, by using QV and UV from one can solve 
for a and /? directly: 



tan{a) 



QV 
UV 



QV 



~UV 



sin(a) cos(q) 



(23) 



(24) 



If one writes as x then uses coefficients a = -^^^ 
and b = ^° one can numerically solve a quadratic for 
and estimate of 7: 



= + l)x^ - {2ab)x + (a^ - 1) 



(25) 



The sign ambiguity in solving for 7=cos^^(x) is re- 
solved by comparing the sign of the computed Mueller 
matrix elements with those of the rotation matrix de- 
rived with these estimated Euler angles. The estimate 
shown here only utilized QQ, QV and UV but provided 
a guess for minimization routines typically accurate to 
better than 2° in our data set. 

In our two-step method, we next solved for the best fit 
Euler angles by doing a normal non-linear least-squares 
minimization by searching the (a, /3, 7) space around our 
initial guesses for minima in the summed squared error: 



estimated starting points but give the same best-fit Eu- 
ler angles to better than one part per million. Note that 
by choosing Euler angle guesses in the proper quadrants, 
the derived rotation matricies avoid the ambiguity of (a, 
13, 7) with ( -180°+a, -13, -180°-|-7). The differences be- 
tween rotation matrix elements and Mueller matrix ele- 
ments for our observations typically much less than 0.1. 
The two-step solution gives us the same Euler angles as 
the direct least squares solution typically better than 0.2° 
and allows for computationally inexpensive and easy im- 
plementation of the least-squares minimization. 




Alpha 



M 100 IM 
BMO 



Alpha - +90 
90 100 ISO 



Fig. 10. — This Figure shows the s qua red error (e^) for the two- 
step solution derived using Equation 1261 Each panel shows a slice 
through the volume (a, (5, 7) at four values of a. Error increases 
from dark to light colors. Good solutions are found at o ~ and 
at a = ±180. 



In order to verify this method, we computed the error 
functions E for both the direct solution and the two-step 
solution and verified the minima with a direct search of 
the entire Euler angle space. We created a grid of Euler 
angles from -180° to 180° (using 2 degree increments) 
and computed the error E as in Equation [TlJ Figure 
shows the derived errors as functions of /3 and 7 for four 
different values of a. Each panel is color coded to show 
the square error E increasing from dark to light. The Eu- 
ler angles found using the two-step method are identical 
to those found with the direct solve within the crude 2° 
sampling we chose for this brute-force search. The error 
functions are always continuous with only two minima in 
the search volume showing convergence is easily achieved 
using any minimization scheme. 



2 3 

E(a, /3, 7) = E E [M'^- - ^)]' (26) 

By using the built-in IDL routines POWELL or 
AMOEBA, we can construct the error function and find 
the minima. Both of these minimization routines require 
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6. APPENDIX: PLANNING SKY OBSERVATIONS USING 
TIME DEPENDENCE FOR MODULATION 

In order to use this technique efhciently, one must con- 
sider the noise propagation when inverting a sequence 
of sky observations to derive telescope Mueller matrix 
properties. There is an analogy between using the time- 
dependent Rayleigh sky to measure polarization proper- 
ties of the a telescope and the retardances chosen to cre- 
ate an efficient modulation scheme for polarization mea- 
surements. The underlying mathematics is the same and 
one can easily derive the requirements on the observa- 
tions needed to derive high-accuracy telescope Mueller 
matrix measurements. Effectively, one is attempting 
to derive properties of a matrix through what different 
groups call demodulations, inversions or deprojections. 
We outline here the analogy between polarimetric modu- 
lation and calculating noise properties with the Rayleigh 
sky to determine a quality observing sequence for deriv- 
ing telescope Mueller matrices. 

Polarimeters modulate the incoming polarization state 
via retardance amplitude and orientation changes. This 
retardance modulation translated in to varying inten- 
sities using an analyzer such as a polarizer, polariz- 
ing beam splitter or crystal blocks such as Wollaston 
prisms or Savart plates. Thes e m odulation schemes 
can vary widely. For example , iCom pain et al. 199 



del Toro Iniesta & Collados"! I200C1| . De Martino et al 
2003, Nagaraju et al. 200'a and pRmczvk et al. 2010 
overview optimal schemes, error propagation and out- 
line schemes to maximize or balance polarimetric effi- 
ciency and create polychromatic systems. There have 
been many implementations of achromatic and polychro- 
matic designs in both stellar and solar commu nities (cf. 
IGisler et ah 1120031 iHmlaoka 112004 IXu et al. ''2006). In 
the notation of these studies, the instrument modulates 
the incoming polarization information in to a series of 
measured intensities (1^) for i independent observations 
via the modulation matrix (Oij) for j input Stokes pa- 
rameters (Sj): 



/ ^ +%/3 +\/3 \ 



O 



1 

t 



1 L L L 



1 

± 

'V3 



1 



V3 I 



(30) 



One recovers the input Stokes vector from a series of 
intensity measurements by inverting the modulation ma- 
trix (O) provided it is a square and multiplying this in- 
verse by the measured intensities. If the matrix is not 
square, one can simply solve the over-specified system of 
equations via the normal least squares formalism: 



O^O 



(31) 



Even for non-square matrices, we can define a demod- 
ulation matrix that captures the transfer properties of 
the modulation scheme: 



= [O^O] 



(32) 



In our case, the Rayleigh sky input parameters become 
the modulation matrix (O^j — Rij) and the formalism 
for noise propagation develop ed in many studies such as 
Idel Toro Iniesta fc ColladosT[2000l apply. 

As an example, we will use the Jan. 27*'* 2010 Mueller 
matrices derived in Figure [SJ The sky polarization ro- 
tated by about 33° during this period and the three cal- 
culated Rayleigh sky inputs scaled to 100% polarization 



are: 



Qr3 Urs 



-0.980 +0.201 
-0.847 +0.534 
-0.716 +0.698 



(33) 



If each measurement has the same noise a and there 
are n total measurements then the noise on each demod- 
ulated parameter (ai) becomes: 



I, = o 



U ^3 



(27) 



This is exactly analogous to our situation where we 
have changed the matrix indices to be i independent 
Stokes parameter measurements for j different sky input 
Stokes parameters: 



2 2 

CT, — na 



Ed 



(34) 



And the efficiency of the observation becomes: 



(28) 



In most night-time polarimeters, instruments choose a 
modulation matrix that separates and measures individ- 
ual parameters of the Stokes vector: 



O, 



1 
1 
1 



\ 

* 

-1 

-1 



0+1 

V 1 -1 / 



(29) 



Other instruments choose only four measurements: the 
minimum number of exposures required to measure the 
Stokes vector. In these schemes, one can balance the 
efficiency of the measurement to minimize the noise on 
each Stokes parameter: 



(35) 



For instance, the normal modulation sequence of equa- 
tion[29]used by most night time spectropolarimeters gives 
n=6 and ef = [l,i,i,-^]. The efficiency balanced scheme 
gives the same relative efficiencies as the normal modu- 
lation scheme but uses only 4 exposures instead of 6. In 
the case of our Jan. 27*'' observations considering only 
qu terms our demodulation matrix is: 



Bij = [R^R]^iR^ 



+1.23 +0.16 
-1.49 +0.43 



-0.48 
-1.54 



(36) 



With this demodulation matrix the efficiency for u 
is only 60% worse than q and we have efficiencies of 
e, = [0.43,0.26] computed with Equation |35l The 
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Mueller matrix derived from these Rayleigh sky obser- 
vations will have similar noise properties between Q and 
U terms. 

One must take care with this technique to build up 
observations over a wide range of solar locations so that 
the inversion is well conditioned. The path of the sun 
throiighoiit the day will create regions of little input sky 
Stokes vector rotation causing a poorly constrained in- 
version with high condition number. For instance, at 
our location in the tropics the sun rises and sets with- 
out changing azimuth until it rises quite high in the sky. 
We are constrained to observing in early morning and 
late evening with the dome walls raised since we may 
not expose the telescope to the sun. This causes input 
vectors at east-west pointings to be mostly q oriented 
with little rotation over many hours. Observations at 
other times of the year or at higher solar elevations are 
required to have a well conditioned inversion. One can 
easily build up the expected sky input polarizations at a 
given observing site with the Rayleigh sky polarization 
equations. Then it is straightforward to determine the 
modulation matrix and noise propagation for a planned 
observing sequence to ensure a well-measured telescope 
matrix with good signal-to-noise. 
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7. APPENDIX: LOVIS SPECTROPOLARIMETER 
CHARACTERIZATION 

In order to expand the capabilities of our spectropo- 
larimeter, we have adapted the optics for use at low spec- 
tral resolution. The polarimeter unit for the spectro- 
graph is a Savart plate providing dual-beam analyzing 
mounted behind the entrance slit with either two rotat- 
ing achromatic retarders or t wo liquid-crystal v ariable re- 
tarders (LCVRs) outline d in Harrington fc K uhn 2008; 
[Harrington et al. 1 120101. We have modified the HiVIS 
echelle housing to allow a flat mirror to be mounted by- 
passing the echelle without disrupting the optical path. 
We call this new mode LoVIS. 




Spectrol Pixel 

Fig. 11. — This Figure shows example LoVIS data when ob- 
serving a spectrophotometric standard star on Sept. 4*'* 2009. 
Wavelength increases from right to left. The Na D lines, Hq and 
the atmospheric A-band can be seen by eye near spectral pixels 
2300, 1400 and 900 respectively. 



With a fiat mirror replacing the echelle, the only dis- 
persive element in the spectrograph is the cross-disperser 
working in 1** order. Only a single order is imaged on 
the detector making the illuminated region of the device 
much smaller. Figure [Til shows a raw data frame from a 
spectrophotometric standard star (HR7940) taken with 
LoVIS in spectropolarimetric mode using the Apogee de- 
tector. There is a single spectral order seen with orthogo- 
nally polarized spectra being displaced roughly 70 pixels 
spatially by the Savart plate. Wavelength increases from 
right to left with the Hq line seen near spectral pixel 1400 
and the 6870A atmospheric band seen near pixel 900. We 
have two separate cross-dispersers availa ble on the rota- 
tion stage blazed for 6050A and 86OOA ( iHodapp et al. I 
[2OOOI iThornton et al. \ \200^. 

The new LoVIS mode gives us much higher sensitivity 
and allows for much faster calibrations at spectral res- 
olutions of 1,000 to 3,000 depending on the slit. The 
1.5" slit had spectral resolutions of 850 and 990 derived 
from Thorium- Argone (ThAr) lines at 5700A and 7000A. 
These lines were spectrally sampled with roughly 9.5 
pixels per Gaussian Full-Width-Half-Max at each wave- 
length. Th e IDL reduction scripts we wrote for HiVIS 
outlined in [Harrington fc Kuhn 1 120081 were adapted to 
this new low-resolution data. Figure [12] shows the ex- 
tracted spectrum from the raw data of Figure [11] 
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Fig. 12. — This Figure shows an example extracted spectrum 
of the spectrophotometric standard star HR7950 taken on Sept. 
4*'' 2009. The H^, Na D lines and the atmospheric A-band are 
identified. 



This change in resolution allows us to observe much 
fainter targets or observe sources at a much higher ca- 
dence. When observing the daytime sky with AEOS and 
LoVIS, we can achieve a spectropolarimetric precision of 
0.1% with a 5 second exposure time using the 1.5" slit. 
The polarimetric images sequences finish in one minute 
and the bulk of the overhead is from rotating the achro- 
matic wave plates in the typical 6-exposure sequence. 

To calibrate this new polarimetric mode and the tele- 
scope polarization we used LoVIS to gather a wide range 
of sky and calibration observations from August 2009 to 
February 2010. All observations were taken with the 
Apogee detector described in Harrington et al. 201(J 
We have sky observations taken solely at the Zenith as 
well as observations on a grid of altitudes and azimuths. 
On several occasions we did observations for 3 or 4 hours 
continuously before sunset to accumulate a wide range 
of polarization inputs at multiple pointings. 

In order to decouple the LoVIS and AEOS polariza- 
tion measurements, we did a number of tests with our 
polarization calibration unit mounted at the slit as well 
as at the entrance port to our instrument at the calibra- 
tion lamp unit. We did tests both without and with the 
image rotator inserted in the the path. This image rota- 
tor is the source of most of the cross-talk in the LoVIS 
instrument. Removing it makes the polarization prop- 
erties of LoV IS much more benign as shown in Figures 
15 and 16 of [Harrington et al. |[20Tol There is a trans- 
missive window that separates the coude rooms from the 
central optical path to the telescope. This window is 
directly after the last fold mirror of the telescope, m7, 
and can be set to either BK7 or infrasil. This window 
can also have a polarimetric effect on the beam and both 
windows were also tested. Effectively we find that the 
polarization properties derived for HiVIS are similar to 
LoVIS and that we can efficiently reproduce polarization 
measurements with both rotating achromatic wave plates 
and LCVR modulators. 

We also measured the induced polarization of the en- 
tire system when illuminated by a point-source utilizing 
unpolarized standard stars. As an example, we show ob- 
servations of many unpolarized standard stars at many 



Mueller Matrices Using Sky Polarization 




5500 6000 6500 7000 7500 

Wovelength In Angstroms 



Fig. 13. — This Figure shows the telescope induced polarization. 
We measured quv spectra for many different unpolarized standard 
stars at many different pointings observed with LoVIS on Sept. S*** 
2009. Each Stokes parameter is shown in percent. The induced 
polarization is at the few percent level with some dependence on 
wavelength. 

pointings from September 5*^ 2009 in Figure [T31 The 
measured polarization is typical of results obtained on 
any individual night. Consistent with our previou s HiVIS 
measurements from [Harrington fc Kuhn 1 120081 the in- 
duced polarization is around a few percent, has some 
mild chromatic properties and does depend on telescope 
focus, the seeing and stellar location on the slit. 

In summary, the polarization properties of LoVIS are 
similar to HiVIS. The induced polarization is below a 
few percent. The polarization calibration optics give 
very pure inputs that are reproduced with both achro- 
matic wave plate and LCVR modulators. This new low- 
resolution mode works efficiently and allows us to obtain 
calibration observations at a much faster cadence. 
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